Generalized unitary Bogoliubov transformation that breaks fermion number parity 
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The standard Bogoliubov transformation is generalized to enable fermion number parity breaking. 
The new transformation can diagonalize fermion Hamiltonians that are quadratic in fermion and 
number parity operators. This new variational freedom cannot lower the energy of Hamiltonians 
that conserve number parity. However, it has value as a minimal mechanism for continuously varying 
between states of different fermion number, regardless of number parity. 

PACS numbers: 73.22.Gk 



The Bogoliubov transformation is a well-established 
mechanism for generating wavefunctions of indefinite 
fermion number. However, these wavefunctions continue 
to have definite number parity. A generalization that can 
break number parity exists but it loses the simple 
structure of a linear transformation between elementary 
fermion operators. Preserving this structure is desirable 
for theoretical considerations of number parity violations 
on equal footing with number conservation violations. 

A new generalization of the Bogoliubov transformation 
is defined by a set of quadratic fermion Hamitonians, 



H = hijdlcj + \9ijc\c\ + \g* l c l c 3 + \Jic\tt + 



7r = exp(i7rnj), (ttci = —cj-k & n = 1) 



(1) 



hu = 



in second quantization notation with hi 
h*jj, and gu = —gji- Lowercase indices are summed 
implicitly over all fermion sites. This model includes a 
fermion number parity operator, tt, with the algebraic 
properties of a Majorana fermion operator but without 
a corresponding site within the Hilbert space. 

As in the Bogoliubov case with // = 0, calculations 
are enabled by an operator set that forms a Lie algebra, 



Q = {1, tccj , c\tt, ctCj, c\cj, c\c'j}. 



(2) 



Any Hcrmitian operator spanned by this set, H — X+X^ 
with X £ span(Q), can be diagonalized by an operator 
of the form: U = exp(Y - Y"f) with Y £ span(Q). The 
transformation applied to the basic fermion operators is 
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(3) 



Z is constrained to have SO(2n + 1) group structure, 



which is a real antisymmetric eigenvalue problem [2| with 



H 



E : 



Im(ff/j + hu) Re(gu - hu) Im(/j) 
Re(/i/j + gu) Im(hu - gu) Re(//) 
-Im(/, 7 ) -Re(fj) 

-ei8u 

ei > 0, 



eiSu 




(5) 



then U diagonalizes H: U^HU = \{ha — £j) + £in,-. The 
transformed vacuum state, U\0), is the ground state. 

The result can be factored into U — U-qUb, where Ub is 
a standard Bogoliubov transformation = 7/ = 0) and 
Ug = exp(#7rc — 9&t\) breaks number parity by coupling 
to a single fermion mode, c := /3iCj, which acts as 
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+ ccos 
(c 1 +c) sin 29 + tt cos 26 



i7rsin26' 



(6a) 
(6b) 



The transformed vacuum that results is a superposition 
of even and odd number parity states, 



U\0) = U B \0) cos9- U-bc ] \0) sm6. 



(7) 



For a Hamiltonian that conserves number parity, there 
is no matrix element between the even and odd states, 
which means that the variational ground state will have 
definite number parity. Multiconfigurational calculations 
where each wavefunction has freedom to independently 
break number parity Q will also have this null outcome. 
There remains the possibility that nonvariational models 
of quantum correlation may artificially stabilize a parity- 
broken ground state. 

Number parity breaking can be a useful theoretical tool 
even when a parity-broken ground state is not stable. For 
example, consider a system where the Hartree-Fock and 
exact ground states have different number parity, 



Hi = \5\hi + A + V(c|4c2 g i 



c[clc 3 c 4 ) 



(8) 
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, det(Q) = l, (4) 



for A = i|A|(l — fi2 + fi3 + hi). If Ue with c = t\ is 
applied to the Hartree-Fock ground state, C2I0), its total 
energy increases to a value of 



for a real orthogonal matrix Q. If Q solves HQ = QE, 



E HF (9) = |<5|sin 2 6>. 



(9) 
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For a post-Hartree-Fock state of the form C/jj^flCjl^) with 
a unitary transformation, U2, that doesn't mix number 
sectors, the variational ground state energy is lowered to 



E vaT (0) = 



W + \A\ 



V 2 



(10) 



The = correlated ground state is unstable to parity 
breaking even though the minimum energy occurs for a 
state of fixed number. This avoids the need to consider 
multiple number-constrained Hartree-Fock solutions as 
reference states for correlation models based on a pure 
state. For correlation models based on a mixed state, 
orbital occupations can be continuously varied without 
breaking number parity. 

In conclusion, the Bogoliubov transformation can be 
generalized to unify the breaking of number conservation 
and number parity. The immediate consequences of this 
generalization are modest, but it may have applications 



to improving the stability of quantum correlation models 
based on a pure reference state. 
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